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ABSTRACT 

Hot core molecules should be detectable in external active galaxies out to 
high redshift. We present here a detailed study of the chemistry of star- forming 
regions under physical conditions that differ significantly from those likely to be 
appropriate in the Milky Way Galaxy. We examine, in particular, the trends in 
molecular abundances as a function of time with respect to changes in the rel- 
evant physical parameters. These parameters include metallicity, dust:gas mass 
ratio, the H2 formation rate, relative initial elemental abundances, the cosmic 
ray ionization rate, and the temperature of hot cores. These trends indicate 
how different tracers provide information on the physical conditions and on evo- 
lutionary age. We identify hot core tracers for several observed galaxies that 
are considered to represent spirals, active galaxies, low-metallicity galaxies, and 
high-redshift galaxies. Even in low-metallicity examples, many potential molecu- 
lar tracers should be present at levels high enough to allow unresolved detection 
of active galaxies at high redshift containing large numbers of hot cores. 

Subject headings: astrochemistry — ISM: molecules — galaxies: starburst — 
galaxies: high-redshift — galaxies: spiral — stars: formation 



Introduction 



Hot cores are commonly observed in regions of massive star formation in the Galaxy, 
through their molecular line emissions. These hot cores correspond to a relatively small, 
dense, warm zone of gas surrounding a newly-formed star. They show a characteristic 
chemistry, distinct from that of molecular clouds in more quiescent regions of the interstellar 



medium (Walmsley & Schilke 1993). Hot cores are rich in relatively large organic molecules. 



some are products of surface chemistry. This characteristic chemistry is generally attributed 
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to a prolonged evolution of the gas and dust at very low temperatures, high densities, and 
large visual extinctions, and it is therefore inferred that the hot cores represent a sample 
of material from the pre-stellar gas that was not incorporated into the newly-formed star 



(Walmsley & Schilke 1993). 



While stellar spectra give information on local elemental abundances after star for- 
mation has occurred, molecular line intensities from hot cores can - with the use of an 
appropriate hot core model - give information on the local elemental abundances before star 
formation occurred. Such an opportunity may be particularly useful in external galaxies 
where the metallicity and other physical parameters may be considerably different from 



those in the Milky Way Galaxy. Hence, Lintott et al. (2005) suggested that hot cores could 



be used as probes of high redshift galaxies in which the rate of star formation and the 
hot core population are much larger than in the Milky way Galaxy. Ideally, one would be 
able to distinguish between different models of stellar evolution for the yields of metallicity 



from zero- metallicity stars. Lintott et al. (2005) showed that if hot cores in a high redshift 



galaxy were sufficiently numerous, then molecular line emission from a large population of 
unresolved hot cores in high redshift galaxies should be detectable. However, there exist at 
present no models of the chemistry of hot cores in the physical conditions likely to pertain 
in high redshift galaxies. 

Therefore, our intention in this work is to calculate the expected chemistry arising in hot 
cores in galaxies which have a wide range of physical conditions that may be substantially 
different to those found in the Milky Way Galaxy. The physical parameters that significantly 
affect the chemistry include gas density and temperature, metallicity, dust : gas ratio, cosmic 
ray flux and its ionization rate, interstellar radiation fleld, etc. Of course, hot cores are only 
one component giving rise to molecular line emission in external galaxies. Other sources of 
molecular line emission are photon dominated regions (PDRs) and cool molecular clouds (see. 



e.g. Garcia-Burillo et al. 2006, 2007). However, hot cores have a characteristic chemistry 
which should give rise to a particular line emission signature. We shall in future work 
examine the emission signature from other types of sources. 



We use as the basis of our computations the hot core model of Viti & Williams ( 1999 ) 



which was developed for the description of hot cores in the Milky Way. In the present work 
we are particularly concerned with galaxies at high redshift in which frequent massive star 
formation and the most intense hot core molecular line emissions occur. The rise time of 
the stellar temperature from that of the cold gas to the Main Sequence is relatively short 
( Bernasconi fc Maeder||199"6 Hanson|[l998 ) and, for present purposes, we shall regard this as 
instantaneous (unlike the treatment of Viti et al.|2004" , in which the evolving chemistry of the 
hot core during the warm-up phase was explored). The model we use describes the physical 
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and chemical evolution of the star-forming region during the long pre-stellar collapse and 
the consequences of the abrupt warming after the star is formed. This two-phase model is 
described in detail in Sect. M 

We select a range of physical parameters intended to cover the range likely to be ap- 
propriate for external galaxies. The metallicity and the dust:gas mass ratio may differ 
substantially from the values of the Milky Way depending on galactic mass and starburst 
activity. The cosmic ray ionization rate is poorly constrained, though possibly related to the 
formation rate of massive stars, and so may be large in active galaxies. Therefore, we must 
choose a wide range for all the relevant parameters; our choice is described in Sect. [3j Hot 
cores in the Milky Way Galaxy are remarkable for the very high visual extinctions (typically 
hundreds of visual magnitudes) associated with them. In galaxies in which the dust:gas mass 
ratio is much smaller than in the Milky Way, the extinction associated with hot cores will be 
smaller. However, a minimum extinction is required for a hot core to exist, so this provides 
a constraint on the choices of density distribution and physical size of a model hot core. 

Our computations describe the evolution of hot core chemistry, and we use them to 
explore the sensitivity of the molecular relative (to the total number of hydrogen atoms) 
abundances to the physical parameters adopted. Our results are presented in Sect. |4j 
Trends of the variation of hot core chemistry with changing physical parameters are identified 
in that Section, and these are the main results of this paper. In Sect. |5] we adopt a set 
of values of physical parameters that may be appropriate for several galaxy types, and we 
predict hot core tracers for these types. Section [6] gives a brief discussion and conclusion. 



2. Model description 



Hot core models have been developed by various authors to describe the evolution of 
chemical abundances of various species within a collapsing cloud during star formation. We 
use the hot core model developed at University College London (hereafter called the UCL hot 
core model; see |Viti fc WiUia"5^[l999| |Viti et al.|[2004| ). To summarize, this model is divided 
into two phases; the first describes the modified free-fall cloud collapse during which time- 
dependent gas phase chemistry occurs and is modified by the gas-grain interaction; it largely 



follows the formulation of Rawlings et al. (1992). The second phase describes the chemical 



evolution after the massive star switched on. The UCL hot core model contains a number 
of free parameters, and our choices are described in detail in Sect. [3j For these massive 
stars, the switch-over from phase 1 to phase 2 of the UCL hot core model is represented by 
an abrupt increase of the temperature (for gas and grains) from 10 K (in the collapse phase) 
to 300 K (in the hot core phase). In phase 2, grains and gas are assumed to be similarly 
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heated and the icy mantles are evaporated. Chemical evolution in the gas phase continues, 
although the dynamical evolution is assumed to have ceased. In both phases of the UCL hot 
core model, the chemical network, is based on more than 1700 chemical reactions taken from 
the UMIST data base ( [Millar et alll997| |Le Teuff et aI]|2000D involving 176 species of which 



42 are surface species. Obviously, all the molecules so far detected in galactic hot cores (e.g. 
Blake et al.|[l987l [Millar fc Hatchell|[l998| [Hatchell et aL][l998[ [Fontani et aL[[2007| have been 



included in this chemistry network together with all the important and familiar molecular 
species that are known from many other studies to play a significant role in the chemistry 
of interstellar gas. The relevant surface reactions included in this model are assumed to be 
only hydrogenation, allowing chemical saturation where this is possible. Note that the hot 
core chemistry necessarily includes species such as CO (not characteristic of hot cores), and 
molecules such as H2O and NH3 which are present in hot cores but undetectable. However, 
these hydrides do become detectable when deuteration to HDO and NH2D has occurred. We 
have not included deuteration (nor any other isotopologues) in the chemistry, in order to 
limit the size of the chemical network. Therefore, in the discussion of results we shall exclude 
water and ammonia, even though it seems clear that these deuterated molecules should be 
detectable. Deuterated species will be considered in later paper. Hydrogenation of carbon 
monoxide on surfaces to methanol is also included. The efficiency we have adopted for this 
reaction is chosen so that, for Milky Way parameters, the observed fraction of CH3OH is 
obtained. The rate coefficients in current chemical network imply that methanol must be 
formed in surface reactions. 

Note that the results reported here are molecular abundances rather than emission 
intensities. The radiative transfer calculations to estimate intensities would require more 
detailed models of hot cores and their populations than is included in our model, and the 
intensities would be significantly model dependent. We do not feel that such computations 
are justified in this exploratory study. 



3. Parameter choices 

The parameter choices to be made in running the hot core model are listed in Table [2[ 
together with the values assigned to these parameters when representing typical hot cores 
in the Milky Way (see Table [T]). We have explored the response of hot chemistry mainly 
to changes in metallicity (with and without consequent changes in the dust:gas mass ratio 
and H2 formation rate coefficient), in cosmic ray ionization rate, in hot core temperature, in 
elemental relative abundances, and - since sulfur-bearing species can be important tracers 
in Milky Way hot cores - in independent variations in the sulfur abundances. We have 
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run over 60 hot core models, but present here resuhs from 19 computations; the parameter 
choices made for these 19 calculations are summarized in Table [2l This table shows that 
we have explored the consequences of reducing solar metallicity by factors of up to 10^, 
and independently of reducing the gas: dust ratio and the H2 formation rate coefficient by 
similar factors. We have also explored the effect of replacing the solar metallicity by values 
predicted by early Universe models, or of reducing the S/H ratio by a factor of 10^. We have 
considered the consequences on hot core chemistry of raising the cosmic ray ionization rate by 
an order of magnitude, and of increasing the hot core temperature from 300 K to 500 K. We 
have considered the consequences of some multiple parameters changes. Where not stated, 
parameters values from Model (representing hot cores in the Milky Way) are used. Note 
that for all 19 computations, rii (initial density), Uf (final density), Ti (cloud temperature) 
and / (UV radiation field intensity) adopt the Model values shown in Table [l] There 
is one proviso to note about these parameter choices. For a hot core to form, the visual 
extinction associated with it must be sufficiently optically thick at far-UV wavelengths, so 
that stellar UV is trapped, warms the core and causes the ices formed during the collapse 
phase to desorb. For Model 13 (which assumes a very low dust:gas mass ratio) it is unclear 
as to whether this criterion is met. If the physical dimension of the hot core is comparable to 
those in the Milky Way (~ 0.03 pc) then the visual extinction for Model 13 would be about 
2.6 mag, equivalent to about 8 mag in the far-UV if the dust grains have similar optical 
properties to those in the Milky Way. 



4. Sensitivity of chemical abundances to variations in the physical and 

chemical parameters 

The main aim of this work is to study the sensitivity of chemical abundances in models 
of massive star formation to variations in physical and chemical parameters which may be 
characteristic of galaxies at different redshifts. A collection of detectable tracers of extra- 



galactic hot cores will be given in Subsect. 4.5 but first we simply analyze the trends that 



the chemistry follows as each parameter was varied. Figures [T] to [8] show examples of 
the most interesting changes while Tables [5} [6} [7] and |8] summarize the temporal trends 
of fractional abundances with variation in metallicity, in cosmic ray ionization rate, and in 
initial elemental abundances. In the following discussion, we rather arbitrarily assume that 
molecules with fractional abundances less than about 10^^^ are probably undetectable, a 
criterion that is roughly satisfied in the Milky Way. 
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4.1. Sensitivity to variations in metallicity 

Models where changes in metaUicity alone (Models to 3 in Table [2]) as well as (more 
realistically) changes in metallicity coupled with other parameters (Models 11 to 13 in Table 
[2]) have been investigated. The first set of models were run in order to isolate the effects due 
solely to metallicity changes (see Figs. [l| g and Tables g g- 

Reducing metallicity leads to no great surprises for most species in that their abundance 
is simply reduced accordingly; however, one notes that as the metallicity is reduced to 1% of 
solar the fractional abundances of some species such as CS tend to decrease with time (with 
the released atomic sulfur going into SO); while for the solar case and 1/lOth of the solar it 
continues to increase with time. 

However, some species do not seem as sensitive to metallicity changes: for example, the 
fractional abundances of CH3CN for different metallicities tend to converge after ~ 10^ yrs, 
making this species rather insensitive to metallicity and therefore (potentially) a good tracer 
of hot cores at high redshift. This 'convergence' is due to a faster destruction (due to cosmic 
rays) at high metallicities (where more ions are present). SO and H2S also tend to converge. 
However, SO conversion into SO2 by cosmic ray-induced reactions is inevitably slowed down 
at low metallicity, so that SO may be more abundant at late times for lower metallicity than 
solar. 

A good discriminator of metallicity is H2CO. Its abundance varies proportionally to 
metallicity as metallicity is reduced from solar values, and then drops below detectable 
values for metallicities less than 1/lOOth of solar. Another interesting species is HCN whose 
fractional abundance is directly proportional to metallicity up to ~ 10^ yrs but then this 
behavior is reversed, i.e. its abundance becomes larger at lower metallicity. This is probably 
due to an increase in ionization for lower metallicity, as can be seen from the behavior of 
HCO"*". The latter, not usually a hot core tracer, is in fact probably only detectable at very 
low (~ 0.1% of solar) metallicity. 

Perhaps more significant may be the trends among models where the metallicity varia- 
tions are also applied directly to the dust:gas mass ratio and the H2 formation rate. However, 
we find that most chemical abundances behave as in the simpler models where only metallic- 
ity varies (see Tables [sjand |6]). It is worth noting that OCS and SO2 abundances are fairly 
insensitive to metallicity changes and are good tracers of hot cores provided the metallicity 
does not drop below about 1% of solar. 
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4.2. Sensitivity to variations in the cosmic ray ionization rate 

The foUowing resuhs concern Models and 14 (defined in Table [2]). Our results are 
shown in Fig [3] and summarized in Table [7} The effects of increasing the cosmic ray 
ionization rate are complex and in general lead to a chemistry that approaches steady-state 
more quickly (see Fig. [s]). For the sulfur-bearing species, a higher ( implies a faster 
destruction of CS and hence a faster increase in the abundances of its daughters species 
such as SO, SO2 and OCS. Nitrogen-bearing species such as HCN and CH3CN increase 
with ( while oxygen bearing species such as CH3OH, H2CO and CH2CO decrease. Some 
species, such as water and CO, are not significantly affected by the larger ionization rate, 
since conversion to these species by both surface and gas-phase chemistries is complete in 
Phase 1. It is interesting to note how oxygen and nitrogen bearing species are anti-correlated 
as a function of (, probably due to the high abundance of ionized carbon in the case with 
higher ionization rate, and this reacts quickly with many oxygen-bearing species. The main 
reservoir of nitrogen is N2 which, when ionized, gives N"*", a reactive species that promotes 
the formation of nitrogen-bearing molecules. 



4.3. Sensitivity to variations in temperature 

We ran two models in which the hot core temperature, T2, was increased from 300 K 
to 500 K (Models 15, with solar metallicity, and Model 16 with 1/lOth of solar metallicity). 
In order to isolate the temperature variations effects, we first compare two solar metallicity 
models where only the temperature is varied (see the top two plots of the Fig. |4]). As 
expected, several species are affected by an increase in hot core temperature: H2S, HCN 
and CH3CN all survive for much longer making them good tracers of temperature. CH2CO 
abundance, on the other hand, is 2 orders of magnitude lower in the high temperature case: 
this is due to its destruction by cosmic ray induced photo-dissociation which has a significant 
temperature dependence. 

At lower metallicity (see the bottom plot of the Fig. |4]), it is interesting to note that 
II2S remains high despite the lower abundance of metals and in fact it is the most abundant 
sulfur-bearing species up until 10^ years when CS becomes comparable: hence, for a higher 
temperature and low metallicity hot core H2S remains a good tracer. 
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4.4. Sensitivity to variations in initial elemental abundances 

The initial elemental abundances (i.e. before star formation occurs) of galaxies at high 
redshifts are by no means well-known. We therefore consider briefly the sensitivity of molec- 
ular abundances in the hot core phase to variations in the initial abundance of the main 
elements (see Models 4 to 10 in the Table [2]). We consider the consequences of adopting 
the elemental abundances arising from stellar evolution models in initially zero metallicity 



gas, as proposed by Chieffi & Limongi (2002); Heger & Woosley (2002); Umeda & Nomoto 



(2002) (see Table [s]). We have also investigated the effect of varying independently the 
initial sulfur abundance. The results for this section are shown in Figs. |5| |6| [7] and |8] and 
summarized in Table [H 

Beside the expected changes in absolute abundances (e.g. H2O abundance increases 
for oxygen-rich models), the most striking variations are, not surprisingly, shown by sulfur- 
bearing species, in particular (see Figs, [s] and [6]) : while for most models CS remains 
abundant for the (canonical) lifetime of the core, in the HW02 model it decreases steeply 
(with atomic sulfur going into SO and SO2) after lO'' yrs. This suggests that in an oxygen- 
rich and sulfur-rich environment CS would not be a good tracer of massive star-forming 
regions while SO, and in particular, SO2 would be enhanced with respect to galactic hot 
cores. The chemical evolution of OCS follows the same trend to late times (~ 10^ yrs) for 
all initial abundances, but its abundance is directly proportional to the amount of sulfur 
initially available so it could be a good indicator of the initial cosmic sulfur abundance. The 
same behavior is found for CH3CN which follows the nitrogen abundance so, if the nitrogen- 
deficient models of Table |3] are correct, this molecule is unlikely to be detectable in high 
redshift hot cores. The sulfur bearing species SO and H2S actually increase in abundance 
after ~ 10^ yrs for Model 6. This is an oxygen-rich model where sulfur is close to solar. 
Finally, although largely insensitive to variations of the initial elemental abundances, H2CS 
is nevertheless always quite abundant and may therefore be used as a tracer of hot cores 
regardless of the chemical composition of the galaxy in question. 

Of other species, H2CO, HON and HNC are of particular interest. The H2CO abundance 
is fairly fiat in its time evolution for oxygen-rich models, at least up to lO^yrs, unless the 
oxygen abundance is very high (as in Model 15). In that case, the H2CO abundance falls 
steadily between lO^yrs and lO^yrs. HON decreases with time for oxygen-rich models, while 
it increases for oxygen-poor models, as HCN is directly destroyed by atomic oxygen. 

In galactic star-forming regions, the sulfur initial elemental abundance is a particularly 
poorly-determined parameter ( Ruffle et al.||1999 Wakelam et al.|2004 ) and yet sulfur-bearing 



species are often used as 'chemical clocks'. The main conclusion of varying solely the initial 
sulfur abundance by 4 orders of magnitude below solar (Models 7 to 10 in the Table [2] 
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the results are shown in Figs. [7] and |8| with summary comments in Table |8]) is that 
non sulfur-bearing species are only sensitive to the initial sulfur abundance if the latter is 
comparable to carbon. So, for example, CH3CN and HCN do not decline as fast when the 
initial abundances of carbon and sulfur are comparable, remaining abundant enough to be 
detected even at late times, while the opposite is true for H2CO which declines after 10^ yrs 
and for CH2CO which never achieves high abundances in Model 0. Sulfur-bearing species are 
of course heavily affected by the initial abundance of sulfur in a fairly predictable way, but 
it is worth noticing that SO and SO2 increase with time for the solar case (Model 0) while 
they decrease in every other model. This may imply that, although commonly considered 
as robust "chemical clocks", sulfur-bearing species may only be useful in this respect for 
star-forming regions closely resembling those in the Milky Way Galaxy. 



4.5. Tracers of extragalactic hot cores 

Tables [5} |6} [7| and [8] show that most of the molecules familiar from studies of hot 
cores in the Milky Way Galaxy are also likely to have significant abundances in hot cores in 
external active galaxies, even though the physical and chemical conditions in those galaxies 
may be very different from those in the Milky Way. We have identified in those tables 
molecules that should trace in hot cores in external galaxies the local metallicity, cosmic ray 
ionization rate, hot core temperature and relative elemental abundances. The molecules are 
CS, OCS, CH3CN, CH3OH, SO, SO2, H2S, H2CS, H2CO, CH2CO, HCN and HNC. Ethanol 
and methyl formate are also known as hot core molecules in the Milky Way Galaxy. However, 
according to our hot core models, they appear to be generally of rather low abundance in 
external galaxies. It is however possible that our models do not account adequately for 
the formation (by surface chemistry) of these two species under conditions appropriate for 
external galaxies. 

For galaxies with metallicities lower than or comparable to that of the Milky Way, the 
molecules with highest abundances (and therefore likely to be useful hot core tracers) include 
CS, SO2, CH2CO and HNC. For metallicities that are at least 1% of that of the Milky Way, 
the following molecules are also be likely to be useful hot core tracers : OCS, H2CS and 
H2CO. The molecule SO can also be abundant, but tends to be an intermediate species in 
the formation of SO2 during the hot core lifetime. We have already noted the convergence of 
the abundance of CH3CN to a value that is rather insensitive to the metallicity; this value 
is rather low (~ 10~^°) but still significant. 

For galaxies with enhanced cosmic ray fluxes but otherwise similar parameters to those 
of the Milky Way, the most abundant species are likely to be OCS, H2CS, CS and SO2; 
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these molecules between them contain most of the available sulfur. Likewise, HCN, HNC 
and CH3CN contain large amounts of the available nitrogen. CH2CO and H2CO have large 
abundances that are enhanced when the cosmic ray flux is raised. 

When the temperature of the hot cores is raised, the ratio of the abundances of HCN 
and HNC approaches unity. 

It is interesting to note that even when the initial elemental abundances are those arising 
from the first generation of stars, hot core molecular abundances can be high. Model 4 shows 
that both CH2CO and H2CO should be prominent tracers. 



5. Selected galaxies and predicted hot core tracers 

In this section, we select several fairly well observed galaxies, representative of different 
morphological types, that may be reasonably compared with some of the model galaxies 
listed in Table [2} O ur aim is to identify the predicted tracers of hot cores in these galaxies. 

We identify Model as a potentially good representative normal spiral. The metallicity 



of IC 342 is estimated to be 1.07 Zg (using 12+log(0/H) = 9.30 from Vila-Costas fc Edmunds 



1992| |Garnett|[l998[ ) while for NGC 4736 we have 1.04 Zq (using 12+log(0/H) = 9.01 from 
Zaritsky et al. 1994). In addition, these two galaxies do not show any bursts of star formation, 
suggesting that a radiation field comparable to that of the Milky Way as well as a hot core 
temperature of 300 K may be appropriate for them. They do not contain any additional 
source of energy such as, for instance, an AGN which could be responsible to an enhancement 
of the cosmic ray ionization rate value. All these physical parameter values suggest that 
Model may be a good representative for these two galaxies. 

A second interesting case is Model 15 having solar metallicity, a standard radiation 
field, hot core temperature of 500 K and normal cosmic ray ionization rate. It could be an 
appropriate model for starburst galaxies such as M 83 (metallicity of 1.06 Zq derived from 



value in Zaritsky et al. 1994) since starburst environments are expected to show higher gas 
temperatures without necessarily any additional sources of energy justifying an enhancement 
of the cosmic ray ionization rate (except in the presence of an AGN which is suggested to 
be a potential source of high energy cosmic rays). Taking a standard value for the radiation 
field for this kind of source may not be appropriate but, as the metallicity is solar, the gas 
to dust mass ratio is 100, the final H2 number density is 1x10^ cm~'^, the visual extinction 
associated with a hot core is ~ 500-600 mag. Therefore, the radiation field has negligible 
influence on hot core tracers in such galaxy types. 
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Galaxies of the type of ICIO are also interesting to study. Garnett (1998) found a value 
1 



of IC 10 metallicity ~ Zq. This source is usually considered as one of the closest starburst 
sources in the Universe even if IC 10 has been recently recognized as only marginally starburst 
( Hidalgo-Gamez|2005 )■ Adding the presence of numerous Wolf-Rayet stars in IC 10, we could 
however suppose a temperature higher than in quiet galaxy. Consequently, a model with a 
z=| Zq, temperature of hot core of 500 K, ^ = 1 (no AGN) appears appropriate. As in the 
previous case, the radiation field is largely irrelevant to the hot core chemistry. We suggest 
ICIO could be represented by Model 17. 

It would be interesting to identify a model for high redshift sources such as the Cloverleaf 
(redshift of Zred. ~ 2.6) or APM 08279 (z^ed. ~ 3.9). For these sources, the physical properties 



are unknown. However, molecular detections already obtained (Gao & Solomon 2004 Wu 



et al.|200"5 Wagg et al.|2005 ) suggest that star formation may be active, even if these galaxies 



are affected by their AGNs. Garci'a-Burillo et al. (2006) have suggested that the unusual 



chemistry they have been observed may be due to a combination of star formation and AGN 
(PDR/XDR) chemistries. In an attempt to describe such sources, we adopt in Model 18 a 
high value of the hot core temperature, a high cosmic ray ionization rate, a high value of the 
radiation field intensity, and a relatively low metallicity (all with respect to the Milky Way). 

Table [9] summarizes the typical molecular tracers of extragalactic hot cores which are 
detectable for these four kinds of sources, where the detectability limit is arbitrarily taken 
to be equivalent to a fractional abundance of lxlO~^^. Here our goal was to give to the 
reader some observational clues derived from this modelling work. The table indicates that 
these should be some differences in the appropriate tracers to use for each galaxy type. The 
overall trends (described in Sect. |4]) of those tracers with physical parameters should enable 
a more detailed description of the galaxy to be made. 



In the two lowest metallicity cases represented in Table [9] (i.e. Models 17 and 18) some 
of the relative abundances of potential tracers are surprisingly large. For example. Model 
17 predicts that in ICIO CS, H2S, HCN and HNC should all have fractional abundances ~ 
10^®, while CH3CN is about one order of magnitude smaller. Model 18 predicts that HCN 
and HNC should have fractional abundances ~ 10~^, CS ~ 10~^, and H2CS and CH3CN 
both ~ 10^^. According to the estimates of Lintott et al. (2005) these abundances are 
sufficiently large that unresolved active galaxies should be detectable in these species even 
at high redshift. Observational tests of prediction made in Table [9] will be described in 



a forthcoming paper. Recent observations of Arp220 (Aalto et al. 2002, 2007) find large 



abundances of some hot core tracers such as HC3N. It is yet not clear whether the chemistry 
is dominated by UV, X-rays or other physical processes. Model 1 or Model 16 may best 
represent this galaxy. 



Gravitational lensing of star-forming regions in high z galaxies may significantly affect 
the molecular line intensities observed. However, relative line intensities should be preserved, 
and the hot core 'signature' should still be discernable. 



6. Conclusion 

We have computed the chemical abundances in hot cores with a wide range of physical 
parameters. In particular, we have discussed the sensitivity of molecular abundances to 
changes in metallicity, to the local cosmic ray ionization rate, to hot core temperature, and 
to the relative elemental abundances that may represent the cosmic composition within 
galaxies early in the evolution of the Universe. In all cases that we have examined there is 
a rich hot core chemistry, even when the metallicity is low. The fractional abundances of 



our predicted tracers are at a level that was shown by Lintott et al. (2005) to provide for 
an active galaxy (containing many hot cores) - even at high redshift - a detectable signal, 
though of course unresolved. 

The chemical network within a hot core is complex, and not all species respond in the 
same way to changes in physical parameters. We find, for example, that some species trace 
linearly any changes in metallicity, while others may respond inversely to such changes. 
Many species are enhanced in abundance if the flux of cosmic rays is increased but other 
species are reduced in abundance. 

We have created models of hot core chemistry that describe crudely several reasonably 
well-observed galaxies. We have predicted hot core tracers for these galaxy types. Detection 
of hot core species should - with the use of models such as those presented here - provide 
detailed descriptions of the physical conditions within active galaxies. 
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Table 1: Standard model parameters (Model 0; see Sect, pi 



Parameter 


Symbol 


Typical Milky 






Way values 


Free fall modifier 


B 


0.1 


Initial number density (phase 1) 


Hi 


300 H cm~3 


Final number density (phase 1 and 2) 




1x10^ H cmr^ 


Temperature (phase 1) 


Ti 


10 K 


Temperature (phase 2) 


T2 


300 K 


External UV radiation intensity 


/ 


1 Habing 


Cosmic ray ionization rate 


c 


1.3x10-1^ s-i 


Visual extinction 


Ay 


600 mag 


Gas: dust ratio 


d 


100 


H2 formation rate coefficient 


R 


1.0x10"^^ X VT cm^s- 


1 


Metallicity 




solar values, see Table 


2 
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Table 2: Input parameters of the UCL hot core models (see Sect. [3] and the text and the 
figures in Sect. |4]). These 19 models have been selected because they are the most relevant 
models for studying the infiuence of the input parameters of the UCL hot core model on the 
predicted relative (to the total number of hydrogen atoms n(H)) abundances (n(X)/n(H)) 
of various species. The abbreviation "ST" represents the standard values listed in Table [T] 
while the abbreviations "CL02, HW02 and UN02" are elemental abundances ratios references 
presented in Table [3j The abbreviation "CB-A, -B, -C, -D and -E" correspond to the values 
listed in Table HI 



Model 


Metalhcity" 


Gas-to-dust 


Ini. Elem. 


S/H 




Temp. 


Comb. 






mass ratio 


Abund. ratios 




(Co) 


hot core (K) 


parameters 





1 


100 


ST 


1.4x10-6 


1 


300 


ST 


1 


1/10 


100 


ST/10 


1.4x10"^ 


1 


300 


ST 


2 


1/100 


100 


ST/100 


1.4x10-^ 


1 


300 


ST 


3 


1/1000 


100 


ST/1000 


1.4x10-9 


1 


300 


ST 


4 


1 


100 


CL02 


7.59x10-6 


1 


300 


ST 


5 


1 


100 


HW02 


1.54x10-^ 


1 


300 


ST 


6 


1 


100 


UN02 


4.06x10-^ 


1 


300 


ST 


7 


1 


100 


ST 


1.4x10"^ 


1 


300 


ST 


8 


1 


100 


ST 


1.4x10"^ 


1 


300 


ST 


9 


1 


100 


ST 


1.4x10"^ 


1 


300 


ST 


10 


1 


100 


ST 


1.4x10"^ 


1 


300 


ST 


11 


1/10 


1000 


ST/10 


1.4x10-^ 


1 


300 


CB-A 


12 


1/100 


1x10^ 


ST/100 


1.4x10-^ 


1 


300 


CB-B 


13 


1/1000 


1x10^ 


ST/1000 


1.4x10-9 


1 


300 


CB-C 


14 


1 


100 


ST 


1.4x10-6 


10 


300 


ST 


15 


1 


100 


ST 


1.4x10-6 


1 


500 


ST 


16 


1/10 


1000 


ST/10 


1.4x10-^ 


1 


500 


CB-A 


17 


1/5 


500 


ST/5 


2.8x10-^ 


1 


500 


CB-D 


18 


1/25 


2500 


ST/25 


5.6x10"^ 


10 


500 


CB-E 



" : Z0 = 1 corresponds to solar values of the elemental abundances ratios; ^ : expressed in 
units of C© = 1.3x10"^^ s-\ 
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Table 3: Initial abundance ratios values used in Table |2] The abbreviations "CL02, HW02 



and UN02" refer to Chieffi fc Limongi (2002), Heger fc Woosley (2002) and Umeda fc Nomoto 



(2002), respectively (see Sect. 3). The standard initial abundance ratios values are from 



Sembach & Savage (1996); Sofia et al. (1997); Meyer et al. (1998); Snow et al. (2002); Knauth 



et al. (2003). We did not include in this table the values of the initial elemental abundance 



S /H ratio since they have been already listed in Table ^ 



ST 



CL02 



HW02 UN02 



C/H 1.4x10"^ 

0/H 3.2x10-^ 

N/H 6.5x10"^ 

He/H 7.5x10-2 

Mg/H 5.1x10-6 



1.4x10-^ 

4.54x10-^ 
5.99x10-11 
7.5x10-2 
1.83x10-5 



1.4x10"^ 1.4x10"^ 

1.53x10-3 1.18x10-3 

1.58x10-^ 3.24x10-^ 

7.5x10-2 7.5x10-2 

1.2x10-4 6.18x10-5 



Table 4: Combinations of parameters coupled with the metallicity and used in the Table. |2] 
(see Sect. |3]). 





ST 


CB-A 


CB-B 


CB-C 


CB-D 


CB-E 




581.1 


59.9 


7.8 


2.6 


117.8 


25.2 


(mag) 














ratio of the number 














densities of grains 


1.0x10-12 


1.0x10-13 


1.0x10-14 


1.0x10-15 


2.0x10-13 


4.0x10-14 


to hydrogen nuclei 














H2 form." 
rate coeff. 


1.73x10-16 


1.73x10-1^ 


1.73x10-1*^ 


1.73x10-19 


6.32x10-1^ 


1.26x10-1* 



: The H2 formation rate coefficient is computed at the beginning of the phase 2 of the UCL 
hot core model and it is expressed in cm3s-i. 
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Table 5: Temporal trends of molecular fractional abundances, computed for different values 
of metallicity. Other model parameters have the standard values listed in Table [1} 



Molecule 


Response to metallicity changes 


CO, H2O 


good tracers of metallicity since flat 


CS 


good tracer of metallicity since grows strongly during phase 2 


ocs 


undetectable at lowest metallicity, otherwise, a good metallicity tracer 


CH3CN 


tends to converge at late times in hot core epoch; 




roughly independent of the metallicity 


CH30H 


good tracer of metallicity, but lost at late times 


so 


converted at late time to SO2, most rapidly for high metallicity cases. 




An inverse tracer of metallicity at late hot core epoch 


S02 


detectable at all metallicities; a linear metallicity tracer 


H2S 


at early times a linear metallicity tracer; lost at late times 




by chemical conversion 


H2CS 


a strong indicator of metallicity, ~ (metallicity)^ 


H2C0 


a linear indicator of metallicity 


CH2C0 


a linear indicator of metallicity 


C2H50H 


undetectable 


HCOOCH3 


undetectable 


HNC 


abundant, but poor metallicity tracer 


HON 


abundant, but poor metallicity tracer; chemical conversion 




at late times 


HCO+ 


good inverse tracer of metallicity, highest for lowest metallicity; 




undetectable for solar metallicity 
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Table 6: Temporal trends of molecular fractional abundances, computed for different values 
of the metallicity coupled with other parameters (see Table 111). 



Molecule Response to changes in metallicity, dust:gas mass ratio, 
and H2 formation rate coefficient 



CO, H2O 

cs 
ocs 

CH3CN 

CH30H 

so 

S02 

H2S 

H2CS 

H2C0 

CH2C0 

C2H50H 

HCOOCH3 

HNC 
HON 
HCO+ 



trends as in Table 5 
trends as in Table 5 
undetectable at lowest metallicity, otherwise, less sensitive 
to metallicity than in Table [s 
fairly independent of metallicity 
traces metallicity; declines steeply after lO^yrs 
trends as in Table |5 

independent of metallicity for higher metallicity values 
trends as in Table [s] 
a strong metallicity indicator 
trends as in Table js] 
follows metallicity 
undetectable 
undetectable 
trends as in Table 5 
trends as in Table |5 
trends as in Table 
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Table 7: Temporal trends of molecular fractional abundances, computed for different values 
of cosmic ray ionization rate {(). Comparison of hot core chemistry for C = 1 ^ind C = 10 
while other parameters have the standard values (see Table 



Molecule Response to higher cosmic ray ionization rate 

CO, H2O unchanged between ( = 1 and ( = 10, insensitive tracer of ( 

CS reduced for ^ = 10 case by ~ x 10, good inverse tracer of ( 

OCS increased for ^ = 10 by ~ x 10, linear tracer of ( 

CH3CN early time (~ 1 xlO^ yrs) reduced for ( = 10, but much increased 

at late time (~ 1 xlO^ yrs), tracer strongly increasing with ( 

CH3OH initially high in both cases, but removed by chemistry 

earlier 1 x 10^ yrs) for ( = 10, good inverse tracer of ( 

SO increased significantly for ( = 10, good tracer for ( 

SO2 increased (by ~ x 10) for ( = 10, linear tracer of ( 

H2S increased significantly for C = 10, good tracer for ( 

H2CS no significant change, insensitive tracer for ( 

H2CO reduced for ^ = 10 after ^ 1 xlO^ yrs, good inverse tracer of C 

CH2CO large in both case, unchanged, insensitive tracer for ( 

C2H5OH undetectable in both cases 

HCOOCH3 barely detectable for ( = I, reduced for ( = 10, but inverse tracer of C 

HNC increased for = 10 by ~ x 10, linear tracer of ( 

HCN increased significantly for = 10, good tracer for ( 

HCO'^ not detectable for ( = I, but increased ~ x 10 for ^ = 10 
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Table 8: Temporal trends of molecular fractional abundances, computed for different values 
of initial elemental abundance ratios (other parameters have the standard values, see Table 

0. 



Molecule 


Response to changes in initial elemental abundances 




UU, rl2U 


unchanged between models (follows total carbon abundance) 






always abundant but removed after ~ 1 x 10^ yrs from hot cores 
in cases with highest total oxygen abundance to form SO 






ivJllvJWo tUtcLl oU.liU.i CLU U.liU.cLliL'C 




CHaCN^ 


follows total nitrogen abundance, a good discriminant between 
solar and early Universe types 




CH3OH 


almost unchanged between models, since the abundance 
mainly follows unchanged carbon 






follows total sulfur abundance 






follows total sulfur abundance 






not very sensitive to the total sulfur abundance, removed at late times 


(^ 1 xlO^ yrs) 




not very sensitive to the total sulfur abundance, removed at late times 


(pa 1 xlO^ yrs) 


H2CO 


inversely sensitive to total oxygen abundance 




CHsCO'^ 


high abundances anti-correlated with total oxygen abundance 




C2H5OH 


negligible but anti-correlated with total oxygen abundance 




HCOOCHa'^ 


low abundances and anti-correlated with total oxygen abundance 




HNC'' 


follows total nitrogen abundance 




HON'' 


follows total nitrogen abundance 




HCO+ 


negligible but anti-correlated with total oxygen abundance 





" : Hot cores tracers CO, H2O follow the total oxygen abundance while CH2CO (and 
HCOOCH3) are anti-correlated with total oxygen abundance; : HCN, HNC and CH3CN are 
all potential discriminants between solar and early Universe elemental abundances through 
the nitrogen content; ^ : typical hot core molecules as SO, SO2 and OCS are correlated with 
total sulfur abundance while H2S and H2CS are rather insensitive. 
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Table 9: Detectable tracers of extragalactic hot core chemistry for four examples of models 
likely to be appropriate for representing the four kinds of galaxy described in the Sect. |5] 
and mentioned in the second line of this table. The limit of detectability has been taken to 
be (n(X)/n(H)) = 1 xlO~^^, as is typical for hot core molecules in the Milky Way. Under 
this limit the species are not detectable (symbol x). Otherwise, they are marked with the 
symbol /. 





Model 


Model 15 


Model 17 


Model 18 


Type of 


Spiral Normal 


Starburst 


low-met allicity 


high redshift 


Galaxy 


IC 342, NGC 4736 


M 83 


IC 10 


Cloverleaf QSO, APM 08279 


CS 


/ 


/ 


/ 


/ 


ocs 


/ 


/ 


/ 


X 


CH3CN 


/ 


/ 


/ 


/ 


CH30H 


/ 


/ 


/ 


X 


so 


/ 


/ 


/ 


/ 


S02 


/ 


/ 


/ 


X 


H2S 


/ 


/ 


/ 


X 


H2CS 


/ 


/ 


/ 


/ 


H2C0 


/ 


/ 


/ 


/ 


CH2C0 


/ 


/ 


X 


X 


C2H50H 


X 


X 


X 


X 


HCOOCH3 


/ 


/ 


X 


X 


HNC 


/ 


/ 


/ 


/ 


HON 


/ 


/ 


/ 


/ 


HCO+ 


X 


X 


X 


/ 
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Fig. 1. — Influence of tlie metallicity alone on the logarithm of the relative (to the total 
number of hydrogen atoms n(H)) abundances (log [n(X)/n(H)]) of various species, with 
respect to the time (in yrs, expressed in logarithm scale as well). We plotted only few 
species which show the most interesting behaviors (see text in Subsect. 4.1). Only models 
from Table [2] are plotted. Model is represented by black lines, Model 1 by dash-dot lines 
while Models 2 and 3 are symbolized by dashed lines and dotted lines, respectively. 



- 24 - 




h23 ■ 



10^ 10^ 10^ 10^ 

time (yrs) 

Fig. 2. — Influence of tlie metallicity alone on the logarithm of the relative (to the total 
number of hydrogen atoms n(H)) abundance of H2S (log [n(H2S)/n(H)]) with respect to the 
time (in yrs, expressed in logarithm scale as well). See caption of Fig. [T| 
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Fig. 3. — Influence of tlie cosmic ray ionization rate alone on the logarithm of the relative (to 
the total number of hydrogen atoms n(H)) abundances (log [n(X)/n(H)]) of various species 
with respect to the time (in yrs, expressed in logarithm scale as well). See the caption of the 
Fig. [T]as well as the text in the Subsect. 4^ We plotted two models whose parameters are 
listed in Table [2] : Model is represented by black lines (top) and Model 14 (^ = 10 x C,q) 
corresponds to the dotted lines (bottom). 
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Fig. 4. — Influence of tlie liot core temperature alone on the logarithm of the relative (to 
the total number of hydrogen atoms n(H)) abundances (log [n(X)/n(H)]) of various species 
with respect to the time (in yrs, expressed in logarithm scale as well). See the caption 
of the Fig. [T] as well as the text in the Subsect. 4^ We plotted three models whose 
parameters are listed in Table [2] : Model is represented by black lines (top two plots) 
while Model 15 corresponds to dotted lines (middle two plots). Model and 15 only differ 
in the value of their hot core temperature. The bottom two plots are for Model 16 (dashed 
lines, low-metallicity T2=500 K, see Table [2|. 
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Fig. 5. — Influence of tlie initial elemental abundances ratios alone on the logarithm of 
the relative (to the total number of hydrogen atoms n(H)) abundances (log [n(X)/n(H)]) of 
various species with respect to the time (in yrs, expressed in logarithm scale as well). See 
the caption of the Fig. [l]as well as the text in the Subsect. 4^ We plotted four models 
whose parameters are listed in Table [2]: Model which is represented by black lines, Model 
4 which is symbolized with dotted lines and the Models 5 and 6 which correspond to dashed 

linps anrl HasVi-Hnt linps rpsr>pptivplv 
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Fig. 6. — Influence of tlie initial elemental abundances ratios alone on the logarithm of 
the relative (to the total number of hydrogen atoms n(H)) abundance of various species (log 
[n(X)/n(H)]) with respect to the time (in yrs, expressed logarithm scale as well). See caption 
of Fig. |5} 
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Fig. 7. — Influence of tfie initial elemental abundance S/H ratio alone on the logarithm of 
the relative (to the total number of hydrogen atoms n(H)) abundances (log [n(X)/n(H)]) of 
various species with respect to the time (in yrs, expressed in logarithm scale as well). See 
the caption of the Fig. [T]as well as the text in the Subsect. |4.4[ We plotted four models 
whose parameters are listed in Table [2] Model which is represented by black lines, Model 
7 which is symbolized with dotted lines and the Models 8, 9 and 10 which correspond to 
small da,shed lines, dash-dot lines a,nd larffe dashed lines, resnectivelv. 



- 30 - 




Fig. 8. — Influence of tlie initial elemental abundance S/H ratio alone on the logarithm of 
the relative (to the total number of hydrogen atoms n(H)) abundance various species with 
respect to the time (in yrs, expressed logarithm scale as well). See caption of Fig. 171 



